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The anomalous scaling of correlation functions in the turbulent statistics of active scalars (like 
temperature in turbulent convection) is understood in terms of an auxiliary passive scalar which 
is advected by the same turbulent velocity field. While the odd-order correlation functions of the 
active and passive fields differ, we propose that the even-order correlation functions are the same to 
leading order (up to a trivial multiplicative factor). The leading correlation functions are statistically 
preserved structures of the passive scalar decaying problem, and therefore universality of the scaling 
exponents of the even-order correlations of the active scalar is demonstrated. 



The old riddle of anomalous scaling in generic turbu- 
lent statistics had been solved recently, albeit in the mod- 
est context of passive scalar and passive vector advection 
[1, 0, E|. Motivated by some recent numerical indications 
we propose in this Letter that under generic condi- 
tions a similar understanding can be extended to active 
scalar (or vector) advection by a generic turbulent veloc- 
ity field. The difference between the two problems is that 
passive fields leave the velocity statistics intact, satisfy- 
ing a purely linear problem into which the velocity field 
is injected as an advection term, while active fields influ- 
ence the statistics of the velocity field. This development 
brings us therefore closer to understanding anomalous 
scaling in the nonlinear velocity problem itself. 

Passive fields are advected in turbulence with typical 
dynamics of the form 



d<f>(r,t) 



+ u(r, t) • V0(r, t) = nV 2 4>{r, t) + f(r, t) . (1) 



at 

Here the field (/>(r,t) can be a passive scalar c(r,t) or 
passive vector v(r,t). k is the diffusivity and f(r,t) is a 
white random force of zero mean with compact support 
in fc-space, acting on the largest scales of the order of 
the outer scale L only. The advecting velocity field is a 
generic turbulent field at high Reynolds number with an 
extended scaling range, typically exhibiting anomalous 
scaling. We use the word "generic" here to distinguish 
from model velocity fields with ^-function temporal corre- 
lations H- "Anomalous scaling" means that multi-point 
correlation functions are homogeneous functions of their 
arguments, with exponents that cannot be guessed from 
dimensional analysis. Thus for example the field (f>(r,t) 
has simultaneous multi-point correlation functions 

F^( ri ,r 2 ,- ■ -,r n ) = W ri ,t)^(r 2 ,t)- ■ -^(r^t)) f , (2) 

where pointed brackets with subscript / refer to averag- 
ing over the statistics of the advecting velocity field and 
of the forcing. Anomalous scaling means that 

FW(Ari,.--,Ar„) = A c -FW(ri,-.-,r n ) , (3) 



with having a non-trivial dependence on n. 

In recent work ^ [| it was clarified why and how pas- 
sive fields exhibit anomalous scaling. The key is to con- 
sider a problem associated with Eq. (|l|) which is the 
decaying problem in which the forcing f(r,t) is put to 
zero. The problem becomes then a linear initial value 
problem, = C(j>, with a formal solution 



(r,t) 



j dr'R{r,r',t)<f>{r',0) , 



(4) 



with the operator R = T + exp ( J Q * dsC(s)j , and T + be- 
ing the time ordering operator. Define next the now time 
dependent correlation functions of the decaying problem: 



C (n) ( 



,i) = (0(ri,t)--.0(r n ,t)) 



(5) 



Pointed brackets without subscript / refer to the decay- 
ing object in which averaging is with respect to realiza- 
tions of the velocity field only. As a result of Eq. (|4|) the 
decaying correlation functions are developed by a propa- 
gator (with r = n,r 2 ,- ■ -,r n ): 

CW(r w --,r n ,t) = JdpP^t) C (n H Pl ,--; Pn ,0) . 

(6) 

In writing this equation we made explicit use of the fact 
that the initial distribution of the passive field 0(r, 0) is 
statistically independent of the advecting velocity field. 
Thus the operator V^ p can be written explicitly 

7>g(i) = {R(r 1 ,pi,t)R(r 2 ,P2,t)---R(r n ,p n ,t)) . (7) 

The key finding J|, || is that the operator pos- 
sesses left eigenfunctions of eigenvalue 1, i.e. there exist 
functions ZW '{r\,r 2 ■ ■ ■ r n ) satisfying 

Z(")(r 1 ,-..,r„)= J dpV™(t)zW( Pl ,.--,p n ) . (8) 
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The functions Z( n > are referred to as "statistically pre- 
served structures" , being invariant to the dynamics, even 
though the operator is strongly time dependent and de- 
caying. How to form from these functions infinitely many 
conserved variables in the decaying problem was shown 
in 0]. The functions Z^ n \r) are homogeneous functions 
of their arguments, with anomalous scaling exponents £ n . 
More importantly, it was shown that the correlation func- 
tions of the forced case, F^ n '(r), have exactly the same 
scaling exponents as Z^(r) H. In the scaling sense 



FW(r) ~ Z (n \r) 



(9) 



This is how anomalous scaling in passive fields is under- 
stood. 

Consider next the dynamics of an active field. For 
concreteness we will think about temperature in turbu- 
lent Boussinesq convection, but the ideas are immediately 
generalized to any other type of active field. The equa- 
tion of motion is formally identical to Eq. (0): 



dT(r,t) 
dt 



-u(r, t) ■ VT(r, t) = K V 2 T(r, t) + f(r, i), (10) 



but now the velocity field is affected by the temperature. 
For an incompressible fluid of unit density M , 



du 
~dt 



u ■ Vu = - Vp + vV 2 u + agTz 



(11) 



Here p, v, a, g and z are the pressure, kinematic viscosity, 
thermal expansion coefficient, acceleration due to grav- 
ity and a unit vector in the upward direction respectively. 
The appearance of T in the equation for u is crucial, and 
changes the scaling exponents of u from Kolmogorov (ap- 
proximately) to Bolgiano (approximately) 0. It makes 
no sense now to consider the decaying problem for T; 
as this field decays, the statistics of the velocity field 
changes, and there is very little to say. On the other 
hand, we can learn a great deal from considering the 
forced solutions. We note that both the passive and the 
active forced equations can be solved in the same way: 

+ J dr' drR{r,r' 7 t-r)f{r',r) . (12) 

with a similar expression for the passive scalar. For t — > 
oo, the first term decays to zero, but the second term 
exists. Next we encounter the difference between the 
active and passive case. Computing the average of T, or 
any other odd moment, we get a finite result: 

(T(r,t))= J dr' J dr(R(r,r',t-r)f(r',r)) . (13) 

We cannot decorrelate the random forcing / from the 
operator R which involves the velocity field. The forc- 
ing / is correlated with T, which is itself correlated, via 



Eq. (|1 ip , with u, and therefore the field T has a first and 
higher odd moments: the probability density function 
(pdf) of T will be asymmetric. For the passive case, the 
equivalent Eq. ( |l3| ) can be simplified by decorrelating / 
from R. As the odd moments of / vanish, the passive 
scalar has zero odd moments, and its pdf is symmetric. 
This difference is fundamental. 

We now argue that all the even moments, and all even 
correlation functions of the active and passive fields may 
be identical to leading scaling order. Consider for exam- 
ple the second order correlation function (n, r 2 ) of 
the active scalar: 

F (2) (n,r 2 ) = J dr'dr" J dr'dr" 

(R(r, r', t - r')/(r', r')R(r, r" , t - r")f(r", r")\U) 

As before, in the case of the active scalar the forces can- 
not be decorrelated from the operators. For the passive 
scalar this can be done with impunity, leading to the 
correct form of the correlation function, i.e. 



(0(ri)0(r 2 )}/ = / dpidp, 

x </(Pl)/(P2)} 



. / drV {2) , fi-r) 

- I ri,r 2 \piP2 y ' 



(15) 



It was shown in || how this expression indeed yields the 
forced correlation function in agreement with Eq. (|5]) . We 
observe however that any correlation function (AB) can 
be written in terms of its Gaussian decomposition and 
the "connected part": {AB) = {A){B) + (AB), where 
A = A — (A) . Identify now A with the product of oper- 
ators R in Eq. ([u]) , and B with the product of forcing, 
and realize that Eq.(|l5|) is the Gaussian decomposition 
contribution to Eq. (M). So the solution ( p^[ ) exists as 
one contribution to (H), the other being what remains 
of the correlation, i.e. the part in which the forcing can- 
not be decorrelated from the operator. A similar phe- 
nomenon occurs for all the even order correlation func- 
tions. But wc know from the theory of the passive scalar 
that the solution ( |l5|) and the equivalent ones for all the 
even-order correlation functions are identified with sta- 
tistically preserved structures whose scaling is anomalous 
H . We propose that the contributions to the even-order 
correlation of the active scalar which are identical to the 
even-order correlations of the passive scalar are leading, 
and are therefore expected to be realized in experiments 
and simulations. The expectation is that the connected 
contributions, that are purely "forced" solutions, are ei- 
ther subleading or have the same exponents. This is a 
crucial conjecture that needs to be tested against exam- 
ples. In H it was shown that this is correct for a par- 
ticular 4th order correlation function of the active and 
passive fields. 

To exemplify the generality and plausibility of these 
ideas we examine in detail a model of active and passive 
scalar for which the statistical object can be computed to 



3 



high accuracy. We consider a variant of the shell model 



at 

~dt 

dCn 

at 




FIG. 1: The scaling exponents for the velocity field (circles) 
and the active scalar field (squares). The solid lines are re- 
spectively 3n/5 and n/5 for the velocity and the active scalar 
fields. Shown in the inset is the pdf of z = u„/ (u^) 1 ^ 2 at shell 
n = 14. 

In this model n stands for index of a shell of k-vector 
k n = k h n , with n = 0, 1, • • • , N - 1. We take h = 2, 
and the parameters used in the simulation are a = 0.01, 
5 = 6 = 6=1, ko = l, k = v = 5x 10~ 4 . The number of 
shells is N = 30. 

Without the coupling to T n , the velocity equation has 

— 1/3 

an inviscid unstable Kolmogorov fixed point, u n ~ k n . 
This is changed by the coupling JtJ, and the system of 
equations for T n and u n exhibits an inviscid unstable 
Bolgiano fixed point, u n ~ k n 3 ^ 5 , T n ~ fc„ . The 
chaotic dynamics renders the statistics of the velocity 
strongly non-Gaussian, cf. inset in Fig[l]. The exponents 
for the active scalar are markedly anomalous, whereas for 
the velocity appear closer to normal, see Fig. |l|. 

To demonstrate the first point of the fundamental dif- 
ference between the active and passive fields, we show 
in Fig. | the pdf's of x = 0n/(</ft) 1/2 w h ere $ is f n 
or C n , for n = 14. One clearly sees the symmetry of 
the pdf of the passive scalar, in contradistinction to the 
asymmetry of the pdf of the active scalar. This is typical 
to all n in the inertial range. This is a demonstration 
of the discussion after Eq.(p"3|). For the passive scalar 



studied in ref. : 

(16) 
(17) 
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FIG. 2: The pdf's of the active (solid) and passive (dashed) 
scalars at shell n = 14. Note that the pdf of the active scalar 
is asymmetric. 



the odd moments vanish, whereas for the active scalar 
they all exist. The situation is altogether different for 
the statistics of even moments. To demonstrate the dif- 
ference we plot in Fig. |^ the (typical) pdf of and 
C\ for n = 9 and 14. In plotting we realize that the 
passive scalar is defined up to a constant, so for the pas- 
sive scalar the pdf is plotted for the rescaled variable 
PCl, where (i = (f%)/(C%) « 0.6327. Note that there 
is only one numerical freedom (3, constant for all n in 
the inertial range. We find very close agreement of all 
the pdf's in the inertial range. The identity of the pdf's 
of and translates automatically to the identity 
(up to a constant [3 m ) of the even-order structure func- 
tions S ,(2m) (fc„) = (4%»), where 4> n = f n or C n . This is 
demonstrated in Fig. [|. We see that the 2nd, 4th and 
6th-order structure functions are barely distinguishable, 
with the same scaling exponents in the inertial range. 
Finally, the identity of the statistics of the squares of the 
passive and active scalars transcends structure functions. 
We consider next multi-point correlation functions, and 
in Fig. H compare the correlation functions (T%T% +5 ) and 
(T%T% +5 T2 n ) to their passive counterparts. The conclu- 
sion is that again the multi-point correlation functions 
are indistinguishable once the passive ones are rescaled 
by (3 q where q is the over-all order of the correlation func- 
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= afc„(w n _iT„_i - hu n T n+ i) + 6fc„(M„T„_i - hu n +iT n+ i) - Kk%T n + fS n ,o 
= ak n (ul_ 1 - hu n u n+1 ) + 6fc„(u„u„_i - hu 2 l+l ) - vk 2 n u n + T n 
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with B, which share however a similar coupling of the 
active scalar to the velocity field. The importance of this 
conjecture warrants considerable further work to uncover 
its provisos and delineate its generality. 




FIG. 4: The even-order structure functions (T 2m ) (circles) 
and (p m C 2m ) (squares), with m = 1, 2 and 3, from top to 
bottom. 



tion. In summary, we offered a discussion and detailed 
numerical evidence to support the conjecture that under 
generic conditions the even-order correlation functions of 
active scalars can be understood via the emerging theory 
of Statistically Preserved Structures of the passive scalar 
counterpart. We have by now two examples, together 




io(T,?T n 2 +5 > (circles) and log lo {0> C%C% +S ) 



FIG. 5: Upper: log 
(squares). Lower: log 10 (T^T^+sTfn) 
logio^ 3 ClCl+sCln) (squares). 
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